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ABSTRACT 

We report improved measurements of temperature anisotropics in the cosmic microwave background 
(CMB) radiation made with the Arcminute Cosmology Bolometer Array Receiver (ACBAR). In this 
paper, we use a new analysis technique and in clude 30% more data from the 2001 and 2002 observing 
seasons than the first release (|Kuo et al.ll2004j ) to derive a new set of band-power measurements with 
significantly smaller uncertainties. The planet-based calibration used previously has been replaced by 
comparing the flux of RCW38 as measured by ACBAR and BOOMERANG to transfer the WMAP- 
based BOOMERANG calibration to ACBAR. The resulting power spectrum is consistent with the 
theoretical predictions for a spatially flat, dark energy dominated ACDM cosmology including the effects 
of gravitational lensing. Despite the exponential damping on small angular scales, the primary CMB 
fluctuations are detected with a signal-to-noise ratio of greater than 4 up to multipoles of i = 2000. 
This increase in the precision of the fine-scale CMB power spectrum leads to only a modest decrease 
in the uncertainties on the parameters of the standard cosmological model. At high angular resolution, 
secondary anisotropics are predicted to be a significant contribution to the measured anisotropy. A 
joint analysis of the ACBAR results at 150 GHz and the CBI results at 30 GHz in the multipole range 
2000 < I < 3000 shows that the power, reported by CBI in excess of the predicted primary anisotropy, 
has a frequency spectrum consistent with the thermal Sunyaev-Zel'dovich effect and inconsistent with 
primary CMB. The results reported here are derived from a subset of the total ACBAR data set; the 
final ACBAR power spectrum at 150 GHz will include 3.7 times more effective integration time and 6.5 
times more sky coverage than is used here. 

Subject headings: cosmic microwave background - 

1. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) radiation on angular scales corresponding to multi- 
pole values of I < 1000 have established a "concordance" 
cosmological model characterized by a negligible spatial 
curvature, 5 % baryonic matter, 2 5% dark matter, and 70% 
dark energy (jSpergel et al.ll2006t ), in good agreement with 
results from other cosmic probes ()Tegmark et all l2004t 
iRiess et alj|2004t iBurles et al.ll2001h . On smaller angular 
scales, the CMB anisotropy is exponentially damped by 
photon diffusion. The damping scale is a measure of the 
angular size of the Silk length at the surface of last scat- 
tering dsHI Il968h iBond & Efstathioul 11984 iHu fc White! 
1997). This portion of the CMB power spectrum, known 
as the "damping tail" , provides a consistency check of the 
cosmological model and independent constraints on several 
cosmological parameters. 

At 150 GHz, anisotropies from the thermal Sunyaev- 
Zel'dovich (SZ) effect are expected to dominate over the 
primary CMB fluctuations for multipoles I > 2500. The 
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— cosmology: observations 

fluctuation power of this "SZ excess" depends sensitively 
on the integrated cluster abundance and cluster gas dis- 
tribution. Accurate measurements of the SZ fluctuation 
amplitude provide an independent measurement of the am- 
plitude of matter perturbations, usually characterized by 

cr 8 . 

The Arcminute Cosmology Bolometer Array Receiver 
(ACBAR) is designed to study both primary and 
secondary CMB an isotropies on small angular scales 
(|Runvan et alj (2003). The fundamental features of 
the angular power spectrum in the range I < 1000 
have been characterized by several ground, balloon, 
and satellite-based experiments. The damping of CMB 
anisotropy power at I > 1000 has b een measured by the 
Cosmic Background I mager (CBI) (jPearson et al.l [2003; 
iReadhead et al.ll2004al), ACBAR and the Very Small Ar- 



ray (VSA) ([Dickinson et al.l l2004h . The excellent agree- 
ment between the observed CMB damping tail power spec- 
trum and the theoretical predictions of the ACDM model 
provide compelling evidence that our interpretation of the 
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CMB is correct. 

The high resolution CMB anisotropy measurements 
reported here extend the low-£ results, such as 
those from the WMAP satellite, to a comoving 
scale of ~ 5 Mpc, and provide strong constraints 
on the sh ape of the primo r dial density perturbation 
spectrum dPeiris et al.l 120031; iMukheriee fe Wanel [2003; 



iBridle et al.ll200& iLeach fcLidd lc 2003). The primordial 
spectrum is close to scale- invariant, consistent with predic- 
tions of slow-roll inflation models. Departures from scale- 
invariance provide valuable clues to the physics of the early 
Universe. 

Measurement s of the high-l CMB anisotropy pow er at 
30 GHz by CBI dMason et al.| [2005 iBond et al.|[200l and 
BIMA (|Dawson et al.ll2006l ) detect power in excess of the 
predictions of the standard cosmological model. If in- 
terpreted as the thermal SZ effect produced by clusters 
of galaxies, this excess power corresponds to a value of 



ysis of optical galaxy clustering l|Tegmark et al 



y anal- 
] [200l 



as that is slightly higher than the value found by anal- 
galaxy 

Efstathiou et al. 2002 . This result hin ges upon the ex- 
cess fluctuation power seen by the CBI, which could be 
the result of contamination by a population of low-flux 
flat-spectrum radio sources. This possibility is being in- 
vestigated with a 30 GHz continuum receiver at the Green 
Bank Telescope (iMason et al.l [20051 ). Alternative origins 
of the excess involve non-standard ingredients in the cos- 
molog ical model, such as, primordial voids and ma gnetic 



fields dGriffiths et aT1l2003l : ISubramaman et alll2003h . The 
frequency spectrum of the thermal SZ effect is distinctly 
different from that of the primordial anisotropy. ACBAR 
observes at a much higher frequency than BIMA or the 
CBI (150 GHz versus 30 GHz), and the addition of the 
ACBAR high-i? measurements can be used to constrain 
the possible origin of any observed excess. 

The fir st AC BAR power spectrum, presented by 
iKuo et all (|2004 ) (K04), was produced from the Lead- 
Main- Trail (LMT) analysis of a subset of the first two years 
of ACBAR observations. Three major improvements are 
made in this work over that first release. First and most 
importantly, the power spectrum is derived from the un- 
differenced temperature maps (rather than "Lead-Main- 
Trail" subtracted maps), significantly reducing the uncer- 
tainties from cosmic variance and instrumental noise. Sec- 
ond, we replace the planet-based calibration with a more 
precise calibration based on the flux of RCW38 as mea - 
sured by ACBAR and BOOMERANG (jCrill et al.ll2003h . 
Third, we include two additional CMB fields observed late 
in the 2002 season which were not used in previous pub- 
lications. These fields increase the total data volume by 
~ 30%. 

This paper is organized as follows. In § [5] we review 
the ACBAR instrument and the CMB observation pro- 
gram. The new analysis algorithm for the un-differenced 
maps is explained in § [31 Section § |H is an overview of the 
calibration using RCW38; the details of cross-calibration 
between BOOMERANG and ACBAR are discussed in Ap- 
pendix A. Systematic tests and foreground contamination 
are discussed in § [5] We present the band-power results in 
§ O including a discussion of the scientific interpretation. 
The ACBAR band powers are combined with the results 
of other experiments to place constraints on the param- 



eters of cosmological models in § [71 The main results of 
this paper are summarized in § [5] 

2. THE INSTRUMENT AND OBSERVATIONS 

The ACBAR instrument was designed to be used with 
the Viper telescope at the South Pole to observe CMB 
temperature anisotropics with an angular resolution of 5'. 
The receiver consists of 16 bolometers, coupled to the 2- 
meter off-axis Gregorian telescope through corrugated feed 
horns. The bolometers are cooled to 240 niK by a 3- 
stage He 3 -He 3 -He 4 sorption refrigerator. The beams from 
the array are swept across the sky at near-constant ele- 
vation by the motion of a flat tertiary mirror. The re- 
ceiver was deployed in December 2000, and CMB data 
were taken during the austral winters of 2001, 2002, 2004, 
and 2005. The bolometer array was reconfigured between 
the 2001 and 2002 observing seasons to double the num- 
ber of 150 GHz detectors; details of the instrument con- 
figuration and perfo rmance in each season are given in 
iRunvan et al.l (|2003l ). while details of the CMB observa- 
tions, data reduction procedures, and beam maps can be 
found in K04. 

The results reported here are derived from the 150 GHz 
data gathered in the 2001 and 2002 Austral winters. These 
data come from observations of four independent CMB 
fields, detailed in Tabic [TJ The power spectrum derived 
from the first two fields, CMB2/CMB4 and CMB5, was 
reported in K04. Since then, we have completed the anal- 
ysis of two additional fields, CMB6 and CMB7, observed 
in July and August of 2002. Each of these four widely- 
separated fields was chosen to include a bright quasar, the 
image of which provides an important check of the tele- 
scope pointing over the entire observation period. 

We derive the noise properties of the data from the 
raw time streams. The photon noise and bolometer noise 
are white, Gaussian distributed, and uncorrelated between 
bolometers. On the other hand, "sky noise" associated 
with atmospheric fluctuations can introduce correlations 
between bolometers. This noise component can be de- 
scribed by the Kolmogorov- Taylor theory, which models 
the turbulence as a screen of frozen fluctuations blown 
through the field of view with a certain wind speed. The 
sky noise properties of the South Pole characteriz e d wit h 
ACBAR are reported in detail bv lBussmann et all (|2005[ ). 
In the CMB power spectrum analysis, we calculate the full 
cross-channel correlations and use them as a data cut. We 
disregard data that show cross-channel correlations (after 
removal of up to a 10th order polynomial from each chop- 
per sweep) higher than 5% of the total noise power. 

3. UN-DIFFERENCED POWER SPECTRUM ANALYSIS 

Following the conventions of the first data release, the 
band-powers q are reported in units of [J.K 2 , and are used 
to parameterize the power spectrum according to 



(1) 



where xb£ are tophat functions; xb£ — 1 for £ G £>, 
and xbi — for I B. The observations in 2001 and 
2002 were carried out in a lead-main-trail (LMT) pattern 
in order to facilitate the removal of slowly varying time- 
dependent chopper synchronous offsets. In the analysis of 
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Table 1 
CMB Fields 



Field 


RA (deg) 


dec (deg) 


Area (deg 2 ) 


Time (hrs.) 


Year 


# of detectors 


CMB2(CMB4) 


73.963 


-46.268 


44(28) 


506(142) 


2001(2002) 


4(8) 


CMB5 


43.372 


-54.698 


31 


1656 


2002 


8 


CMB6 


32.693 


-50.983 


29 


351 


2002 


8 


CMB7 


338.805 


-48.600 


32 


420 


2002 


8 



Note. — The central quasar coordinates and size of each CMB field observed by ACBAR during 2001 and 2002. The fifth column 
gives the detector integration time for each field after cuts. This represents approximately 50% of the total time spent observing 
CMB fields. The last column gives the number of 150 GHz detectors. 



K04, we differenced the CMB maps according to the for- 
mula M — (L + T)/2, and derived the band-powers from 
the differenced maps. While this conservative strategy 
can prevent potential systematic errors, it increases the 
effect of instrumental noise and reduces sky information. 
At high-£, where the uncertainties in the band-powers are 
dominated by the instrument noise, the anisotropy power 
signal-to-noise ratio (S/N) derived from a LMT differenced 
map is a factor of 



VLMT 



SN L 



AIT 



SN 



Vl 2 + (1/2)2 + (1/2)2 
y\/2 2 + (V4/2) 2 + (V4/2) 2 



smaller than that for an un-differenced map with an iden- 
tical amount of observing time (ignoring the small corre- 
lations between differenced fields). The higher noise in L 
and T reflects the reduced (50%) observing time on these 
fields compared to the main field. 

ft is possible to generalize the noise-weighted co-added 
map analysis outlined in K04 to avoid this S/N degrada- 
tion, while ensuring that the power spectrum is not con- 
taminated by the chopper synchronous signals. The key 
to this generalization lies in the fact that the band-powers 
can be derived from arbitrary linear combinations of the 
time stream. The linear combinations can be constructed 
such that the un-differenced maps are cleaned of chopper 
synchronous offsets. The results presented in § [5] are de- 
rived from such "cleaned", un-diffcrcnccd maps, using the 
method outlined below. 

Suppose d a (a — l..n t ) are n% time-ordered measure- 
ments of CMB temperature. This vector is the sum of the 
noise component n a and the signal component s a . The 
data from each chopper sweep are filtered to remove chop- 
per synchronous offsets before they are co-added. For ex- 
ample, since the beams approximately follow a great circle 
on the celestial sphere, changing the chopper angle causes 
the detectors to look through a different atmosphere thick- 
ness and produces signals that are functions of the chop- 
per position. Formally, the filtering is done by operating 
on the time stream d with the "corrupted mode projec- 
tion" matrix II. This operation results in the cleaned time 
stream d = lid. 



Large angular scale chopper offsets are largely removed 
by this mode projection. In the analysis of K04, any resid- 
ual small angular scale offsets were eliminated through 
LMT differencing. Alternatively, chopper synchronous off- 
sets can be removed by subtracting the average chopper 
function, determin ed with the w hole map, from the indi- 
vidual data strips (Runyan 2002). However, this assumes 
that the offsets do not change in time and elevation. In 
this paper, we remove a chopper synchronous offset from 
each data strip where the amplitude of the offset at each 
sample in the strip is free to vary quadratically with ele- 
vation in the map. The mean of these quadratic functions 
in elevation gives the mean chopper function; at zeroth or- 
der, the mean chopper function is removed. The quadratic 
variation allows for slow changes in the synchronous offset 
as a function of time and elevation. Projecting out these 
corrupted modes eliminates any detectable residual chop- 
per offset with minimal loss of signal in the final power 
spectrum. The loss of information at high-£ is negligible, 
since the data contain many more degrees of freedom than 
the removed modes do. 

Mathematically, the corrupted mode projection matrix 
n is now the product of two matrices, n = Il^ni. The op- 
erator ni is the original n matrix used by K04 that adap- 
tively removes polynomial modes in RA, where the order 
of the polynomial removed typically depends on the atmo- 
spheric conditions. The additional operator II2 removes 
quadratic modes in DEC independently for each of the 
lead, main, and trail fields. We perform systematic tests 
in S}5]to confirm that any residual offset is smaller than the 
noise level of the final power spectrum. The resulting time 
stream d is then co-added into a map T according to the 
telescope pointing model. Next, a signal-to- noise eigen- 
modc truncation is applied, and only the modes expected 
to have non-negligible signal-to-noise ratio are retained. 
This significantly reduces the computational requirements 
of the analysis. Since both the projection of corrupted 
modes and signal-to-noise eigenmode truncation are lin- 
ear operations, the resulting data vector A is also a linear 
combination of the original time stream, represented by 

A = Ld. 

The noise covariance matrix is given by 

Cat = L(nn*)L*. 
The signal component in the time stream s is the con- 
volution of sky map X(r) with B a (r), the beam function 
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during measurement a, 

s a = J d 2 r1(r)B a {r). 

The signal component of A is 

K 9 = E L » S « = / ^(r)S(r), (2) 

where the pixel-beam function Fi is given by 

F i (v)=J2 L ic t B a [r). (3) 

a 

In the flat sky limit, the theory covariance matrix is given 

by 



C T{ij} = (AiAj) 



d 2 rd 2 r'Fi (r)Fj (r') (T(r)T(r')) 



;C e -e 



il'(r— r ) 



d 2 i 
W) 2 



Cf ^(1)^(1), 



(4) 



where (1) is the Fourier transform of Fi (r) . To perform 
the iterative quadratic band-power estimation procedure, 
it is necessary to know the partial derivative of Ct with 
respect to each of the band-powers q B , which according to 
equation (fi} is given by 

dC T{ij} f d 2 l 2ir Xm p ni) p M (5) 



dq B J (2Tr) 2 e{e + i) 

Note that this algorithm does not assume that the instru- 
ment bea ms stay constan t durin g the observations. As de- 
scribed in iRunvan et alJ (|2003h . the ACBAR beam sizes 
are weak functions of the chopper position. K04 adopted 
a semi-analytic expansion to correct for these effects to 
first order. To verify that the effects due to non-uniform 
beams are small, we developed two end-to-end pipelines. 
In the first pipeline, the pixel-beam functions Fi(r) are 
calculated explicitly during the co-adding process. The 
bandpowers in Table [3] are analyzed with this algorithm. 
In the second pipeline, an averaged beam is used for the 
entire map. The difference in the power spectra from the 
two pipelines is negligible. 

In the analysis of K04, we assumed that the noise is sta- 
tionary in chopper position after LMT subtraction. In the 
current treatment, we relax this assumption and calculate 
the full two dimensional correlation matrix directly from 
the time stream data without using Fourier transforms. 

All the numerical calculations are performed on the 
National Energy Research Scientific Computing Center 
(NERSC) IBM SP RS/6000. The evaluation of Fi(r) and 
its Fourier transform are the most computationally de- 
manding steps in this analysis. After Ct, Cjv and Ct,b 
are calculated, standard likelihood maximizing procedures 
are used to find the band-powers qs and uncertainties 
(Bon d et al.l Fl998l . The results of this analysis are pre- 
sented in Tableland Figure [T] 

4. CALIBRATION 

RCW38 is a compact HII region in the Galactic plane 
at a declination similar to the ACBAR CMB fields. It 
has a large and stable flux and serves as the primary cal- 
ibrator for the ACBAR observations. We determine the 



absolute flux of RCW3 8 using map s from the 2003 flight of 
BOOMERANG (jMasi et alJ pOOl . hereafter B03), which 
are calibrated relative to the WMAP experiment with an 
absolute uncertainty of 1.8%. RCW38 does not have a 
black-body spectrum, requiring spectral corrections for the 
calibration of CMB anisotropics. However the similarity 
in the spectral responses of the 150 GHz bands in the B03 
and ACBAR experiments ensures these corrections to be 
small. Here we outline the calibration procedure, leaving 
the details to Appendix [A"l 

ACBAR typically observed RCW38 before and after 
each CMB observation. Comparisons between the B03 
and ACBAR maps of RCW38 are used to determine the 
absolute calibration of the CMB fields to an uncertainty 
of 6.0%. For roughly 50% of the 2002 season, we observed 
RCW38 with only half the 150 GHz detectors (4 out of 
8). During these periods, the RCW38 calibration was ap- 
plied to the remaining detectors by comparing CMB power 
spectra derived from each half of the detectors. The cali- 
bration of the CMB4 field (observed in 2002) is extended 
to the 70% overlapping CMB 2 field (observed in 2001) 
by comparing power spectra from each field. In the first 
ACBAR release, the 2001 and 2002 data sets were cali- 
brated with an accuracy of 10% using observations of Mars 
and Venus respectively. We determine the corrections to 
this planet-based calibration to be 0.911±0.072 for CMB2 
(2001) and 1.128±0.066 for CMB4-7 (2002) in CMB tem- 
perature. The 2002 observations dominate the final power 
spectra, and the final results have essentially the same 6% 
temperature calibration uncertainty as the 2002 data. 

5. SYSTEMATIC UNCERTAINTIES AND FOREGROUNDS 

5.1. Jackknife Tests 

We performed a series of tests to constrain the ampli- 
tude of potential systematic errors in the power spectrum 
results. As described by K04, the "first half minus second 
half jackknife is a very powerful test for time dependent 
errors, such as a changing calibration, inconsistency in the 
beam or pointing reconstruction, and time varying side- 
lobe pickup. In addition, high-<? jackknife band-powers 
constrain the mis-estimation of noise. We perform this 
test on the joint CMB power spectrum and find the band- 
powers of the chronologically differenced maps are consis- 
tent with zero (Fig. 

Similarly, the data can be divided in two halves ac- 
cording to the direction of the chopper motion. Micro- 
phonic vibrations due to the chopper turn-arounds, er- 
roneous transfer function corrections, or effects of wind 
direction could produce a nonvanishing signal in the jack- 
knife band-powers. We find that the power spectrum of 
the left-right differenced maps is also consistent with zero. 

To ensure that any residual chopper synchronous offset 
is below the noise level, we have developed a new system- 
atic test that enhances the contribution of any such offsets 
relative to the CMB. In this test, the band-powers are de- 
rived from an LMT sum map, L + M + T, in which the 
residual synchronous offsets (the same in each field) are 
enhanced relative to the (random) CMB fluctuations by a 
factor of 3 in power (neglecting the small correlations at 
low-£) . The resulting band-powers are compared with the 
un-differenced band-powers to check for systematic devi- 
ations. This test is particularly sensitive to any resid- 
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ual chopper offsets. We find no significant deviation in 
the LMT sum band-powers from the un-differenced band- 
powers. When compared with the model ACDM power 
spectrum, we do notice a slight rise in the LMT sum band- 
powers for £ > 2300. ft is difficult to assess the significance 
or the origin of this low level trend. However, even if it is 
caused by a residual systematic effect, the contribution to 
the joint band-powers would be smaller than the statisti- 
cal uncertainty in the reported band powers in this paper 
after accounting for the factor of 3 amplification. 

5.2. Foregrounds 

At frequencies below the peak of the CMB intensity 
(~ 200 GHz) , the contribution of extra-galactic radio point 
sources to the observed CMB temperature anisotropy de- 
creases rapidly with increasing observing frequency. In 
addition to the negative spectral indices of the majority of 
the radio sources, the flux-to-temperature conversion fac- 
tor, (dB v I dTc m b) 1 > reaches a minimum as the observing 
frequency approaches the peak of the CMB. In particular, 
this factor is nearly 15 times smaller at 150 GHz than at 
30 GHz. The measurements of ACBAR are therefore much 
less susceptible to contamination by radio point sources 
than experiments operating at 30 GHz such as CBI, BIMA, 
and VSA. We construct templates using the positions of 
the known radio source s from the 4.85 GHz Parkes-MIT- 
NRAO (PMN) survey (jWright et al.lll994D . and project 
out their contributions to the power spectrum estimations. 
Using the method described by K04, we remove them from 
the data without making assumptions about their fluxes. 
Of 200 PMN sources in the observed CMB fields, we de- 
tected the guiding quasars and six additional sources with 
significance greater than > 2.8a. These sources tend to 
have shallow, and in some cases inverted spectral indices. 
Table [2] lists the parameters of the PMN sources that are 
detected in the ACBAR fields; the detection threshold of 
> 2.8cr, corresponds to a false detection rate of 1. The un- 
certainties are calculated using Monte-Carlo simulations, 
and are dominated by contributions from CMB primary 
anisotropies. With the exception of the guiding quasar in 
each of the CMB fields, the effect of removing the point 
sources on the band-powers is not significant. 

Thermal emission from interstellar dust also has the 
potential to contaminate the measured power spectrum. 
The ACBAR CMB fields are located in the region s of low 
Galactic dust emission. iFinkbeiner et all (|l999f ) (FSD) 
combined observations from IRAS, COBE/DIRBE, and 
COBE/FIRAS to generate a multi-component dust model 
that predicts the thermal emission at CMB frequencies 
with an angular resolution of 6 arcminutes. We apply the 
ACBAR filtering to the predicted dust maps for 150 GHz, 
and find the expected RMS to be at the /iK level. Assum- 
ing the ACBAR maps contain the FSD dust templates 
with amplitudes parametrized by the quantity £, the ob- 
served maps T can be written as the sum Tqmb +£,Tfsd- 
After cross-correlating the dust template maps Tfsd with 
the observed maps, we find that the Galactic dust is un- 
detectable in the ACBAR 150 GHz data. The 1-a upper 
limit on the amplitude parameter is £ < 2.6, consistent 
with the FSD predictions (i.e., £ = 1). As in the case for 

11 http:/ /www. roe. ac.uk/ifa/shades/ 

12 http:/ /chilel. physics. upenn.edu/blastpublic/ 



the radio source flux measurements, the uncertainty in £ 
is dominated by the CMB primary anisotropies. There- 
fore, dust with the morphology of the FSD maps does not 
significantly contribute to the observed anisotropy power. 
Nonetheless, the reported band-powers are calculated with 
the dust template mode projected out in each of the fields. 

Dust emission from high redshift star forming galax- 
ies can be a significant foreground contaminant in mil- 
limeter waveleng th CMB maps. Recent observations with 
SCUBA/ JCMT (ISmail et all 120021: IChapman et al.ll2002t 
iBorvs et all 120031 ), Bol ocam/CSO (lLaurent et al J (20051 ). 
and MAMBO /TRAM (jGreve et all 120041) . provide con- 
straints on both the source counts and the spectral de- 
pendence of these proto-galaxies. However, despite the 
tremendous progress made in studying these sources, their 
contributions at 150 GHz are still highly uncertain. The 
uncertainties come from the low number statistics in 
source counts and spectral dependence, difficulties in mod- 
eling the survey bias and completeness, and poorly stud- 
ied angular correlations. In the absence of more decisive 
measurements, such as might be produced by the ongoing 
SHADESZ3 survey and the BLAST^ experiment, we ig- 
nore the clustering noise component and estimate the Pois- 
son contribution from these proto-galaxies at the ACBAR 
observing frequencies. 

The SCUBA results constrain the 850/im source counts 
with an uncertainty of ~ 40%. There are considerable un- 
certainties in extrapolating this result to lower frequencies 
because the process depends not only on the dust prop- 
erties, but also on the cosmic star formation history and 
the source evolution. Observations carried out at two dif- 
ferent wavelengths, 1.2 mm with MAMBO, and 1.1mm 
with Bolocam, can pot entially proyide th is extrapola- 
tion phenomenologically. iGreve et all (|2004h find that the 
MAMBO and SCUBA source counts agree if the MAMBO 
counts are scaled up in flux by a factor of 2.5, correspond- 
ing to a spectral dependence of S v oc z/ 2 65 . We scale 
the fi t for the SCUBA source count results (|Borvs et alj 
2003() to 150 GHz (2 mm) accordi ng this spectral depen- 
dence. Using the formulae given in IScott fe Whitel (119991) . 
we find the contribution of these sources to the CMB 
power spectrum at 150 GHz to be V e ~ 37(£/2500)VK 2 . 
Since most of this fluctuation power comes from sources 
with fluxes between 0.1 and lOmJy, the result is insen- 
sitive to a flux cut-off greater than lOmJy. The Bolo- 
cam team ^Laurent et al.l 120051 ) found fewer sources at 
1.1mm than MAMBO did at 1.2 mm, implying a steeper 
spectrum for the sources. A fluctuation analysis of the 
same data also su ggested lower source counts, especially 
fo r fluxes < lmJv (iMalonev et al.| [2005). From Figure 15 
of lLaurent et al.l (|2005f ). we estimate the spectra between 
850 nm and 1.1 mm to scale as S v oc v . Compared with 
the 37/iK 2 at I = 2500 from the MAMBO/SCUBA ex- 
trapolation, which is just below the instrumental noise of 
ACBAR, the Bolocam/SCUBA model results in a negligi- 
ble value of 4/xK 2 . Future sub-millimeter and millimeter 
observations are needed to fully characterize the properties 
of these sources. In the interpretation of the ACBAR data, 
we cautiously assume that the high-€ band-powers are not 
significantly contaminated by high redshift proto-galaxies. 
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Table 2 

Millimeter Bright PMN Sources 



Source Name/Position 


Field 


S4.85 (mjy) 


5150 (mjy) 


a 150/4.85 


PMN J0455-4616* 


CMB2 


1653 


2898 ± 60 


0.15 


PMN J0451-4653 


CMB2 


541 


360 ± 58 


-0.13 


PMN J0439-4522 


CMB2 


634 


383 ± 73 


-0.16 


PMN J0514-4554 


CMB2 


422 


197 ± 68 


-0.23 


PMN J0515-4556 


CMB2 


990 


524 ± 65 


-0.20 


PMN J0253-5441* 


CMB5 


1193 


1799 ± 66 


0.11 


PMN J0210-5101* 


CMB6 


3198 


1268 ± 86 


-0.28 


PMN J0214-5054 


CMB6 


61 


186 ± 63 


0.31 


PMN J2235-4835* 


CMB7 


1104 


656 ± 63 


-0.16 



Note. — These sources from the PMN 4.85 GHz catalog a re detected at > 2.8 a significance with ACBAR, corresponding to a 
false detection rate of 1. The fluxes at 4.85 GHz (5*4.85, from lWright et al.1 HHl)) and 150 GHz (S150, measured by ACBAR) are 
given. The spectral index a is defined as S v oc v a . The uncertainties associated with S150 are dominated by the CMB fluctuations. 
The central guiding quasars (one in each of the 4 fields) are marked with asterisks (*). These sources, as well as the undetected 
PMN sources, are projected out from the data using the methods described by K04 and do not contribute to the power spectrum 
measurements reported in this paper. 



6. RESULTS AND DISCUSSIONS 

6.1. Power Spectrum 

Applying the analysis method described in the previous 
sections to the ACBAR 2001 & 2002 150 GHz data leads 
to the power spectrum shown in Figure [T] A comparison 
with other recent CMB results is shown in Figure EB The 
model curves in both figures are the "WMAP3+ACBAR" 
best fit model. We report the decorrelated band-powers, 
since the description of their statistical properties requires 
fewer parameters. Plotting the decorrelated band-powers 
also simplifies the visual comparison between models and 
the measurements. T he decorre l ation transformations are 
defined according to iTegmarkl (|1997t ) . The same trans- 
formations are applied to the window functions, which 
conve rt a model Ct to the theoretical band-powers (jKnoxl 
1999). Following K04, we use the offset lognormal func- 
tions (jBond et al.ll2000l ) to fit the likelihood functions, and 
report the fit parameters q, a, x for each band. The band- 
powers, uncertainties, and lognormal offsets are given in 
Table [3J this information along with the corresponding 
window functions are available for download from the 
ACBAR website 0. 

The ACBAR data are consistent with the results of 
other CMB experiments, and fit the model predictions 
for a flat, A-dominated universe with a low baryon den- 
sity. A narrow peak is clearly seen in the power spec- 
trum at £ ~ 820, corresponding to the third harmonic of 
the acoustic oscillations in the early universe. This de- 
tection is in agreement with previous detectio ns of this 
feature by BOOMERANG ([Jones et all [20061 ) and fur- 
ther c onfirms the coherent origin of the cosmic perturba- 
tions (|Albrecht et al.| [l996). Despite the nearly exponen- 
tial damping at the high-£, the primary CMB fluctuations 
are detected with a signal-to-noise ratio of greater than 4 
up to £ = 2000. The photon diffusion mechanism predicted 
at high-l 7 is verified to a high degree of accuracy. 

13 http:/ /cosmology.berkelcy.edu/group/swlh/acbar/index.html 



6.2. Anisotropics at £ > 2000 

Various theoretical models predict that the secondary 
anisotropy induced by the thermal SZ effect in clusters of 
galaxies begins to dominate over primary ani sotropy at £ > 
2000 for standard cosmologica l parameters (jCoorav et al.l 
I200fi iKomatsu k Seliakll2002f ). The level of the signal is 
extremely sensitive to the normalization of the the matter 
power spectrum, usually parameterized by the present-day 
RMS mass fluctuation on 8 hr 1 Mpc scales, a&. The SZ 
effect has a clear spectral signature. In the nonrelativistic 
limit, the thermodynamic temperature difference from SZ 
effect (ATsz) is given by 

AT SZ ( e* + l \ 

= y I x — — 7 - 4 > (6) 



where x = 



, . x 

Tcmb V e x - 1 

h » = z,/56.8GHz. 



kT c 



The quantity y is 



known as the Compton parameter and is proportional 
to the integrated electro n density along the line of sight 
through the cluster (e.g., iPeacockl (|l999j i). The CBI Deep 
field observations at 30 GHz detect power on scales cor- 
responding to I > 2000 in excess of the predicted pri- 
mary CMB anisotropy. This "excess power" has been 
interpreted as the SZ effect produced by intervening 
galaxy clusters dMason et al.ll2003t iReadhead et alj|2004al 
iBond et al.ll2005l ). On the other hand, a variety of models 
including non-standard primordial e ffects have also been 
proposed as possible exp lanations (|Grifhths et al.l 120031 : 
ISubramanian et al.l [20031) . The unique photon emission 
spectrum of the thermal SZ effect distinguishes it from 
these alternative explanations for the observed anisotropy. 

The ACBAR band powers corresponding to the small- 
est angular scales lie slightly above the best fit WMAP3 
ACDM model. The four highest £ bins jointly produce an 
excess of 51 ±42/xK 2 after the model primary power spec- 
trum is subtracted. The combination of this result with 
measurements at lower frequencies can be used to con- 
strain the photon emission spectrum of the excess, shed- 
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ding light on its origin. We perform a joint analysis of 
the CBI results and the new ACBAR data at £ > 2000, 
assuming the contributions from primary anisotropy are 
known. In each experiment, the theoretical band-powers 
for primary anisotropy are calculated from the product 
of the ACDM power spectrum and band window func- 
tions, which are then subtracted from the observed band- 
powers. A two-dimensional likelihood function is calcu- 
lated from these excess band-powers and their uncertain- 
ties, where the two parameters are the ratio of the 30 GHz 
and 150 GHz excess, £, and the power at 30 GHz, a 2 
(in iiKq MB ). We then marginalize over the cr| param- 
eter and plot the likelihood function for the power ratio 
£ in Figure [U Since ACBAR measures significantly less 
power at 150 GHz, the data disfavor sources that result 
in a blackbody spectrum (i.e., ( — 1). Using the ACBAR 
and CBI frequency response and equation ([5]) , we calculate 
the power ratio ( = 4.3 for the thermal SZ effect. From 
the likelihood plot, we conclude that it is 4.5 times more 
likely that the excess seen by CBI and ACBAR is the re- 
sult of the thermal SZ effect = 4.3) than a primordial 
process (C = 1). Since the expected ratio of flux at 30 GHz 
to 150 GHz from radio sources is expected to be <0.1, such 
sources are disfavored as being responsible for the excess 
power seen by CBI at about the same significance (~ 1.2a) 
with which ACBAR detects excess power. Additional data 
from ACBAR or other higher frequency instruments will 
be required to make a definitive statement about the origin 
of the excess power seen by CBI and BIMA. 

7. COSMOLOGICAL PARAMETERS 

In this section, we estimate cosmological parameters for 
a minimal inflation-based, spatially-flat, tilted, gravita- 
tionally lensed, ACDM model characterized by six param- 
eters, and then investigate models including extra param- 
eters to test extensions of the theory. For our base model, 
the six parameters are: the physical density of baryonic 
and dark matter, fl^h 2 and Vl c h 2 ; a constant spectral index 
n s and amplitude In A s of the primordial power spectrum, 
the optical depth to last scattering, r; and the ratio of 
the sound horizon at last scattering to the angular diame- 
ter distance, 9. The primordial comoving scalar curvature 
power spectrum is expressed as V s (k) — A s (k/k n )^- ns ~ 1 \ 
where the normalization (pivot-point) wavenumber is cho- 
sen to be k n — 0.05 Mpc -1 . The parameter 6 maps an- 
gles observed at our location to comoving spatial scales 
at recombination; changing 9 shifts the entire acoustic 
peak/ valley and damping pattern of the CMB power spec- 
tra. Additional parameters are derived from the basic set. 
These include: the energy density of a cosmological con- 
stant in units of the critical density, Q\; the age of uni- 
verse; the energy density of non-relativistic matter, f2 m ; 
the rms (linear) matter fluctuation in S/i^Mpc spheres, 
as ; the redshift to reionization, z re ] and the value of the 
present day Hubble constant, Hq, in units of kms _1 Mpc~ 1 . 
Tilted primordial spectra indicate the presence of a tensor- 
induced anisotropy component, however, we do not include 
this potential contribution due to its uncertain amplitude. 
The influence of the tensor component would only be sig- 
nificant at low-^, not in the regime which ACBAR probes. 
We also restrict this work to flat ACDM models, motivated 

14 http: / / cosmologist.info / cosmomc 



by the observed curvature being so close to zero. However, 
we have run models with non-zero curvature fix , and find 
that they reproduce the standard geometrical degeneracy 
associated with Vtx and 

We have also considered two extensions to the basic 
model which could potentially impact the interpretation 
of the ACBAR bandpowers. These extended models in- 
clude flat ACDM models with a running scalar spectral in- 
dex characterized by the derivative dn s /d In fc(fc„), and flat 
ACDM models with a Sunyaev-Zel'dovich contribution to 
the angular power spectrum with amplitude parametrized 
by a . We also investigate a model where both a run- 
ning spectral index and a SZ contribution are considered 
simultaneously. 

The parameter constraints are obtained using a Monte 
Carlo Markov Chain (MCMC) sampling of the multi- 
dimensional likelihood as a function of model parame- 
ters. Our software i s based on the pub licly available 
CosmoMCQ package (jLewis fc Bridlell2002j) . CMB angu- 
lar power spectra and m atter power spectra are computed 
using the CAMB code (jLewis et al.lboOOD . We approx- 
imate the full non-Gaussian ba ndpower likelihood s with 
an offset lognormal distribution (jBond et al.l l2000h found 
by explicit fits (see K04 for a detailed discussion of the 
calculation). Our standard CosmoMC results inclu de the 
effects of weak gravitat ional lensing on the CMB (|Seliakl 
1996| iLewis et al.ll200d[) . Lensing effects in the tempera- 
ture spectrum are expected to become significant at scales 
£ > 1000, hence it is important to include this effect when 
interpreting the ACBAR results. The major effect of lens- 
ing is a scale-dependent smoothing of the angular power 
spectrum which diminishes the peaks and valleys of the 
spectrum. Inclusion of lensing in the model improves the 
fit to the data for all experiment combinations. However, 
we find that the parameter mean values and uncertainties 
are largely unaffected by the inclusion of lensing with some 
exceptions, in particular the introduction of lensing tends 
to increase the value of <7g. 

The typical computation consists of 8 separate chains, 
each having different initial random parameter choices. 
The chains are run until the largest eigenvalue of the 
Gelman- Rubin test is smaller than 0.1 after accounting 
for burn-in. Uniform priors with very broad distribu- 
tions are assumed for the basic parameters. The standard 
run also includes a weak prior on the Hubble constant 
(45 < H < 90 km s _1 Mpc -1 ) and on the age of the 
universe (> 10 Gyrs). We also investigate the influence of 
adding Large Scale Structure (LSS) data from the 2 degree 
Field Galaxy Redshift Survey (2dFGRS) dCole et al.ll2005D 
and t he Sloan Digital Sky Survey fSDSS) (jTegmark et akl 
2004). When including the LSS data, we use only the 
band powers for length scales larger than k ~ 0.1/iMpc -1 
to avoid non-linear clustering and scale-dependent galaxy 
biasing effects. We marginalize over a parameter b 2 which 
describes the (linear) biasing of the galaxy-galaxy power 
spectrum for L+ galaxies relative to the underlying mass 
density power spectrum. We adopt a Gaussian prior on b 2 
centered around b g = 1.0 with a conservative width equiv- 
alent to Sb g — 0.3; all parameters except as and r are 
insensitive to this width. 
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7.1. Base Parameters Results 

The results for the basic flat tilted ACDM parame- 
ters are shown in Table 01 The confidence limits are 
obtained by marginalizing the multi-dimensional likeli- 
hoods down to one dimension. The median value is ob- 
tained by finding the 50% integral of the resulting like- 
lihood function while the lower and upper error limits 
are obtained by finding the 16% and 84% integrals re- 
spectively. The CMBall data combination includes: the 
ACBAR results presented here; the WMAP 3 year TT, 
TE, an d EE spectra, with th e EE not included at higher 
I as m iHinshaw etall (120061): the CBI extended mosaic 



2.3 J 

100/6» = 95.9±£;£, 



results dReadhead et al.l 2004af) and polar ization results 
(|Readhead et alj|2004bt [Sievers et alJl2005D. combined in 



the manner described in S ievers etaUljSooa P 5 ! the DASI 
two year results (lHalverson et al.|2002f) : the DASI EE and 
TE bandpowers ftLeitch et al.ll2005l ): the VSA final results 
([Dickinson et al.l 120041) ; the MAXIMA 1998 flight results 
(jHananv et al.ll2000D : and the TT, T E, and EE results 
from the BOOMERANG 2003 flight (iJones et alj [2001 
iPiacentini et alj|2006t iMontrov et al.ll2006[) . Only £ > 350 
bandpowers are included for BOOMERANG because of 
overlap with WMAP3 (although inclusion of the lower £ re- 
sults leaves the parameter results essentially unchanged). 
While ACBAR and BOOMERANG are both calibrated 
through WMAP, this is a small contribution to the total 
uncertainty in the ACBAR calibration and we treat the 
calibration uncertainties as independent in our parameter 
analysis. Although the DASI, CBI and BOOMERANG 
2003 EE and TE results for high I polarization are in- 
cluded, they have little impact on the values of the pa- 
rameters we obtain. 

In all our runs we have used the u pdated WMAP3 likeli- 
hood code (http://lambda.gsfc.nasa.g ov/) which includes 
an up dated point-source correction cf. Huffe nberger et ail 
(2006]) and foreground marginalization on large angular 
scales. These updates result in a small incre ase in the 17 m 
and a s values compared to those reported in lSpergel et al.l 
(|2006h . 

The results for the basic model parameter set with var- 
ious combinations of data are summarized in Fig. [5J The 
most striking feature of the results is that the solutions de- 
termined from WMAP3 alone are quite compatible with 
the extension by ACBAR (and that of the other data) to 
higher £. This consistency means that the additional CMB 
data (including ACBAR) have little impact on the cos- 
mological parameters determined by WMAP 3. We have 
tested the effect of a significantly smaller ACBAR cali- 
bration error, such as we anticipate for the final ACBAR 
release. We find a much larger impact on the parameter 
values and errors; the values are similar to those found for 
CMBall+LSS. 

With the original ISpergel et all (|2006h WMAP3 likeli- 
hood code, there was a shift in as and Q m to higher val- 
ues when additional data was included. However, with 
the updated WMAP3 likelihood code, the addition of the 
ACBAR and CMBall bandpowers leads to essentially no 
shift in erg and fi TO ; however, including the LSS data does 
still result in a slight increase in these parameters. The 
new likelihood code corrects the lower power in the third 



acoustic peak which was leading to low values for erg and 

The comoving damping scale, determined as a de- 
rived cosmological parameter using only the ACBAR and 
WMAP3 data is R D = 10.5 ± 0.2 Mpc. The correspond- 
ing angular scale is Id — 1355T. 5 . These values for Rjj and 
£d are in excelle nt agreement with values obtained using 
earlier datasets ()Bond et al.l 1200 3). We also find the co- 
moving sound crossing distance is R s = 147.8l 2 3 Mpc 
with a corresponding angular scale > 
in agreement with the value for 9 in Table [4j 

Inclusion of lensing in our standard parameter runs in- 
creases the best-fit model likelihoods in all cases. The 
difference between the log likelihoods of the lensed and 
non-lensed models for the WMAP3 run is Ami = 0.86. 
The log likelihood difference increases to 1.7 with ACBAR 
included, 2.46 with CMBall, and 3.69 for the CMBall+LSS 
data combination. The mean values of the parameters do 
not shift significantly with the inclusion of lensing; for ex- 
ample, a 8 increases from 0.778 to 0.788 for the CMBall 
data set and from 0.804 to 0.813 for CMBall+LSS. The 
best-fit T>gs for the lens and no-lens cases look quite sim- 
ilar, but the subtle smoothing of the peaks and troughs 
by lensing results in a better fit to the the data for each 
combination of experiments. 



7.2. Running Spectral Index 

The first release of WMAP data showed evidence for 
running of the CMB power spectrum spectral index, 
articularly when c ombined with measurements of LSS 
Spergel et al.l l2~003). Extending the basic model to allow 
for running of the spectral index around the pivot point 
K = 0.05MPC" 1 yields dn s /d\nk(k±) = -0.053^fg for 
WMAP3 only. The tendency for negative running in- 
dices is due mostly to the low I end, where the multi- 
poles are lower than the standard ACDM model. The 
contribution from the high I end is less significant. Since 
the WMAP3 results extend to reasonably high £, the ad- 
dition of the ACBAR results shifts the constraints only 
marginally dn s /d\nk(K) = -0.045 j^g- The effect 
of adding the ACBAR data can be seen most clearly 
in Fig. [7] which shows the correlation between n s and 
dn s / 'dlnfc. The central value is similar, but the errors 
are further reduced with the CMBall + LSS combina- 
tion, dn s /dln k(k*) = — 0.047^o 021- Similar to the results 
from WMAP1 and earlier versions of the CMBall data set 
( Bo nd et aill2003l : iMacTavish et al.ll2006f ) , a negative run- 
ning is still favored at about the 2-<7 level by each of the 
data combinations considered. The models including run- 
ning favor significantly lower values of the scalar spectral 
index, n s — 0.903^ ; 028 . However this result depends on 
the choice of pivot point fc*: a smaller value would yield a 
higher result while a higher one would give an even lower 
result. 

7.3. Sunyaev-Zel'dovich template extension 

As described in Section 16.21 fluctuations from the ther- 
mal Sunyaev-Zel'dovich (SZ) effect are expected to dom- 
inate over the damped primordial contributions to the 
CMB anisotropy at multipoles beyond £ ~ 2500. The 

15 We exclude the band powers below t = 600 from the CBI extended mosaic results to reduce the correlation with the TT band powers of the 
CBI polarization dataset which influence the sample-dominated end of the spectrum. 
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magnitude of the SZ signal depends strongly on the over- 
all matter fluctuation amplitude, erg. We have modified 
our parameter fitting pipeline to allow for extra frequency 
dependent contributions to the CMB power spectrum and 
have implemented it in a simple analysis using a fixed tem- 
plate Cf z for the shape of the thermal SZ power spectrum. 
The template was obtained from large hydrodynamical 
simulations of a scale-invariant (n s = 1) ACDM mode l 
with a 8 = 0.9 and Sl b h = 0.029. fSee iBond etail (|2005t ) 
for a detailed description of the simulations.) Recently the 
WMAP team have used a different SZ temp late based on 
analy tic estimations of the power spectrum (|Spergel et alj 
2006). It is characterized by a slower rise in £ than the 
simulation-based one, which cut nearby clusters out of the 
power spectrum. There has been no fine-tuning of the 
spectra to agree with all of the X-ray and other cluster 
data. This may have an effect on shape, especially at high 
L 



We add an SZ contribution Cf z 



sz 



(a^) 2 fXi to 



the base six parameter model spectrum. Here f v 
is the frequency-dependent SZ pre-factor and a sz = 
(Tg^ 2 (fib/i/0.029) is a scaling factor determined from hy- 
drodynamical simulations. We consider two cases: (1) the 

scaling parameter a sz is slaved to a^ 2 (Oft/i/ '0.029); (2) 
a sz is allowed to float freely. Including this SZ template 
with all parameters varying is complementary to the anal- 
ysis of § 16.21 which directly compared the residual CBI and 
ACBAR bandpowers in the excess region. In that more re- 
strictive analysis, the primary power spectrum is fixed and 
/„ is allowed to vary as well as a broad-band excess power. 
We found that the spectrum of the excess was compatible 
with the SZ effect, motivating the SZ-restricted study con- 
sidered here. 

Regardless of the data combination, we find that in- 
cluding an SZ component in the model has little effect on 
the determination of the basic cosmological parameters. 
This can be seen by comparing Table [5] with Table [H We 
also find that it has little effect, whether a sz is related to 
cosmic parameters through a sz = a^ 2 (Qbh/ 0.029) or is 
allowed to float freely. Note that the SZ results break the 
A s e~ 2T near-degeneracy (as does weak lensing, though not 
as strongly.) 

We begin with the combination of the ACBAR and 
WMAP3 data. When a sz is allowed to float freely, we ob- 



tain a 



sz _ 



0.57 



+0.20 
-0.56- 



We can use the above definition to 



map the floating SZ amplitude parameter a sz to a corre- 



0.84t°;23- Th e low significance of excess 



sponding crg SZ ^ 

power in the ACBAR data results in weak constraints on 
fjf 2 , particularly for the lower limit. We have not tabu- 
lated the results for the slaved a case since it results 
in extremely small changes in erg. For example, fits to 
the ACBAR+WMAP3 band powers give a 8 = 0.771^ 
when we include the SZ contribution in the model, and 
eg = 0.781q'q5 when we ignore it. 

When the high I bandpowers of CBI and BIMA are in- 
cluded in the analysis, there is a significant detection of 
excess power. Both the CBI and BIMA bandpowers are 
from 30 GHz intcrfcrometric observations and have higher 
/„ values than ACBAR. For the slaved case, the errors 

16 We note that the non-Gaussian nature of the SZ signal is included 
The effect of the sample variance tends to open up the allowed range 



tighten slightly while the central value remains station- 
ary with er 8 = 0.79lg 04- F° r the floating case we find 

a sz = 0.79l8:g3 which maps to erf Z) = 0.92!°;^. 

The CMBall + BIMA combination results in uncertain- 
ties for (Jg SZ ' ) which are comparable to those of erg. A visual 

summary of the results is shown in Fig. [5] where we plot 
(sz) 

both og and er g against the spectral index for a num- 
ber of data combinations. It is interesting to note that 

the tension between as and erg is relaxed by the inclu- 
sion of the LSS data which increases the value of as- We 
caution, however, that the fit depends on the SZ template 
shape and its extension into the higher I regime probed 
by BIMA. This analysis assumes that no additional fore- 
ground sources, such as dusty proto-galaxies, contribute 
significantly to the observed anisotropy. It is worth 
mentioning that the baseline lensed model for the CMB 
continues to be a better fit to the data regardless of the 
inclusion of the SZ extension. The decrease in log likeli- 
hood for the fit to the model including SZ, when lensing 
is taken into account, is A In L = 2.06 for CMBall+BIMA 
and A In L = 2.38 for CMBall+BIMA+LSS. We also find 
that neglecting lensing in the model tends to increase the 
existing tension between the derived as and a§ z values. 

We find that including the SZ template and running 
spectral index simultaneously selects a more negative run- 
ning index due to the greater downturn in high £ power. 
For example, with the CMBall + BIMA + LSS combi- 



nation we get dn s /d\nk 



-0.60 



+0.022 
-0.02 



1 compared with 



the no-SZ result dn s /d\nk = — 0.471q] u2 i- In either case, 

the excess power leads to crg SZ ' = 0.921^06) virtually un- 
changed from the results for the model without running. 

8. CONCLUSIONS 

We have measured the CMB angular power spectrum up 
to multipole values of £ ~ 3000 using the complete data 
set from the 2001 and 2002 ACBAR 150 GHz observations. 
The data are analyzed with a refined method in which the 
band-powers are calculated from un-differenced maps. We 
calibrate the data by comparing the flux of the galactic 
HII region RCW38 in the overlapped region of ACBAR 
and BOOMERANG; the flux scale of BOOMERANG is 
calibrated through comparison with WMAP. The new cal- 
ibration is found to be consistent with the previous planet- 
based calibration, but with uncertainty reduced from 10% 
to 6.0% in temperature. 

We have carried out various jackknife tests to show that 
the ACBAR band-powers are not compromised by sys- 
tematic errors. In deriving the power spectrum, we have 
projected out known foreground modes including the FSD 
dust template and the PMN radio sources, however, these 
projections have no significant effect on the final band- 
powers. The contribution to temperature fluctuations at 
150 GHz from high redshift dusty protogalaxies remains 
uncertain. However, extrapolating from recent observa- 
tions near 270 GHz, we determine that these sources are 
unlikely to contribute significantly to the ACBAR band- 
powers. 

The band-powers presented in Tableware the most sen- 
sitive measurements of the CMB temperature anisotropy 

in the BIMA results but was not taken into account in the CBI analys is, 
towards lower <rg values (Goldstein ct al. 2003; Rcadhcad et al. 2004a). 
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to date in the range of 1000 < £ < 3000. The power spec- 
trum continues to support a spatially flat ACDM cosmol- 
ogy, with a low baryonic density. Since the WMAP3 data 
now extends into the third peak, the addition of the damp- 
ing tail data from ACBAR results in only minor changes in 
the values and uncertainties of the standard cosmological 
parameters. For all combinations of data we have con- 
sidered beyond WMAP3 alone, the baseline lensed CMB 
model results in significantly better fits than models ne- 
glecting lensing. 

The ACBAR data has also been used to place interesting 
constraints on secondary anisotropies. The band-powers 
for I > 2000 are significantly smaller than those reported 
by the CBI and BIMA experiments and provide only a 
suggestion (1.2<r) of excess power above what is expected 
from primary CMB anisotropy (§ 16. 2[) . The excess power 
has a frequency spectrum consistent with the thermal SZ 
effect and inconsistent with thermal sources such as pri- 
mary anisotropy. Because of the weak detection of excess 
power by ACBAR, radio sources are slightly disfavored as 
the source of the signal, but cannot be ruled out. 

Theoretical work suggests that the thermal SZ effect 
should be the dominant source of secondary anisotropy. 
The expected amplitude of the thermal SZ effect is ex- 
tremely sensitive to as- Adding the SZ amplitude to our 

cosmological parameter runs, we infer values for cr s that 
are somewhat higher than the as found from the standard 
parameter runs, but consistent within the uncertainties. 



This tension is further reduced when LSS data is included 
in the parameter runs. 

The results presented here are derived from a subset of 
the total ACBAR data set which is currently being ana- 
lyzed. The final ACBAR power spectrum at 150 GHz will 
include 3.7 times more effective integration time, 6.7 times 
more sky coverage, and a direct calibration from compar- 
ison with WMAP3. This will result in significantly re- 
duced uncertainties across the entire power spectrum and 
improved constraints on standard cosmological parameters 
as well as secondary anisotropies. 
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APPENDIX 

A. CALIBRATION 



The calibrat ion used for the first ACBAR data release was obtained by observations of Venus and Mars, as detailed in 
iRunvanl (|2002f h In this section, we describe a new calibration based on the comparison of B03 and ACBAR observations 
of the galactic HII region RCW38. B03 was calibrated by an q; m -based c omparison of CMB structure with WMAP in 
the region of B03 sky coverage with a 1.8% uncertainty (|Jones etaT] [20061 ). The B03- ACBAR cross-calibration method 
is described below, with a detailed accounting of uncertainty in Table IA61 



BOOMERANG- ACBAR Cross Calibration ACBAR made daily high signal-to-noise maps of the galactic source 
RCW38. B03 also mapped portions of the galactic plane including RCW38 (Fig. IA.0.0.0[) . allowing a direct comparison of 
the high signal-to-noise maps made by the two experiments. The experiments have different scan patterns, beam widths, 
and spatial filters that can effect the measured flux. We resample the B03 map using pointing information for each ACBAR 
observation to generate an ACBAR-equivalent B03 observation. The ACBAR maps are smoothed to simulate the effect 
of Boomerang's larger beam. Large spatial modes in both experiments are corrupted; in ACBAR by chopper-synchronous 
offsets and in B03 by the high-pass filter. We simultaneously fit a quadratic offset and Gaussian source model from each 
scan of the ACBAR and B03 maps which removes these modes without affecting the amplitude of a point source. After 
coadding the channel maps, we integrate the flux within a 18' radius of the source. The integrated flux is robust to small 
misestimates or changes in the beam size. The measured flux ratio and the associated uncertainties are listed under Ratio 
of BOS over ACBAR in Table [AS 

We use Monte Carlo techniques to estimate the transfer function of this method. Using a model of RCW38 and its 
surroundings, we generate simulated timestreams for the observations with each experiment. Maps are created from the 
timestreams and are filtered as described above. The ratio of the transfer functions is found to be ACBAR/ BQ3 = 
1.056 ± 0.002. We have tested the dependence of the transfer function on the assumed signal template and include a 
3% uncertainty in our calibration due to this effect. This technique is readily adapted to include the effect of the beam 
uncertainty for each experiment and we find that the beam contributes 1.35% to our estimated uncertainty. The effect of 
the transfer function and the associated uncertainty are listed under Transfer Function in Table [A"6l 

RCW38 has a much different spectrum than the CMB, and the effective CMB temperature difference it produces depends 
on the photon- frequency. The calibration described here is based on observations with ACBAR's 150 GHz channels and 
Boomerang's 145 GHz channels which have similar bandpasses (Fig. lA.O.O.Op . We account for the small difference in 
bandpass by conv olving the measured spectral response of each experiment with a model of RCW38's spectrum from 
(M asi et al.[ [2006). If two maps nominally calibrated in CMB temperature units are integrated about RCW38, the true 
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Table 3 

Joint Likelihood Band-powers 



i range 


l eff 


q ( M K 2 ) 


a (/iK 2 ) 


x (/iK 2 ) 


351-550 


428 


2680 


284 


-698 


551-650 


605 


2225 


280 


-218 


651-750 


700 


2018 


227 


-286 


751-850 


804 


2796 


276 


-509 


851-950 


910 


1662 


174 


-257 


951-1050 


1003 


1282 


132 


-111 


1051-1150 


1102 


1284 


124 


-146 


1151-1250 


1204 


1116 


108 


-40 


1251-1350 


1303 


877 


92 


-43 


1351-1450 


1403 


782 


89 


63 


1451-1550 


1502 


563 


73 


118 


1551-1650 


1601 


524 


70 


139 


1651-1750 


1703 


351 


62 


266 


1751-1875 


1810 


254 


54 


272 


1876-2025 


1943 


294 


57 


307 


2026-2175 


2096 


278 


74 


520 


2176-2325 


2242 


59 


68 


508 


2326-2500 


2395 


196 


100 


851 


2501-3000 


2607 


190 


120 


1625 



Note. — Band multipolc range and weighted value £ e ff, dccorrelated band-powers qs, uncertainty <tb , and log-normal offset xb from the 
joint likelihood analysis of CMB2, CMB5, CMB6 and CMB7. The PMN radio point source and IRAS dust foreground templates have been 
projected out in this analysis. 



Table 4 

Basic 6 Parameter Constraints 



libit 

n c h 2 



T 

ris 



log[10 10 A s 
Age [Gyrs 

z r e 
H 



WMAP3 
0226+ 00008 

U.UZZD_q 0007 

n m«+ - 008 
u.iU8_ 009 

n 097+ - 012 
u.u»/_ 014 

-+0.017 
-0.017 
(-0.08 

-•-o.oe 
7^+ - 03 

13.7«' 2 
25+ ' 03 

(-0.06 
-0.05 
}+2.7 
5 -2.6 
7+3.0 
-2.3 



0.9667 
3.051 



0.78: 
11.81 
72.71 



0.092 
0.964 



WMAP3 + ACBAR 
2 2 5+ 0007 

U.UZZO_ 00 007 

ins 4 " - 008 

U.1UO_ 007 

1 042+ 004 

-,+0.014 
0.014 
+0.016 
-0.015 

3.051^7 
0.76lg;» 
13.7«- 2 
0-241H 4 
0.78t°;™j 

- + 2.3 

-2.6 



CMBall 



11.51 
72.9 -;;1 



0.0226 



-0.0006 
-0.0006 
+0.006 



0.110_ .006 

l-042±g;gg| 

0.092ir i4 
0.9631°;™ 

o or+0.06 
J- uo _0.06 

0.75lg;S 
13.71°:! 

n 9^+0 03 
u - z,J -0.03 
n 7O+0-04 
u - ' y -0.04 

n.6l 2 : 2 

72.3 + ?1 



CMBall+LSS 

n226 +0 0006 

U.UZZO_ooo06 

n i-ir+0.005 

U.11Q_ 00B 

1.0421°:°" 
0.090ini 3 3 
0-960±S-°« 

o n7 +0.05 
°- u '-0.06 

^,+0.03 
-0.03 

13.71S-.1 

n oo+0. 03 
U - 28 -0.03 
1-0.03 
-0.03 
-+2.2 
-"-2.4 

70.4l 2 j 



0.721 



0.811 
11.51 



Note. — Results for the basic parameter set. The runs all assumed flat cosmologies, uniform and broad priors on each of the basic six 
parameters and a weak prior on the Hubble constant (45 < Ho < 90 km s _1 Mpc -1 ) and the age (> 10 Gyr). All runs included the effect of 
weak gravitational lensing on the spectra which is significant for the t > 1000 scales probed by ACBAR. The results are nearly identical if we 
include an SZ template, with the largest, but still relatively minor, impact on u% and r. 



calibration factor K will depend on the measured flux ratio 1 303/ Iacbar, the bandpass of each experiment t v , the 
spectrum of RCW38 S* CW38 and the known blackbody spectrum of the CMB %-\t cmb : 

Ibo3 d 
K = * R 

1 ACBAR 

J t^BAn X 2 S HCW 3 S d , J tBOS^^^ 
Jt^S^W3S dp J^CBAR^^^- 
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Table 5 

SZ Template Parameter Constraints 



n r h 2 



T 



log[10 10 A s 
Age [Gyrs 
as 



,T SZ 

°8 



H 



ACBAR+WMAP3 

0224+ 0008 
u.uzzy:_ 0006 

n i06 +0 008 

i-04il8:Si 

0.0961^11 

U -°'-0.56 

962+ ' 018 
u.»uz_ 01(i 

76+ 03 

u - '"-0.04 
iq 7+0.2 
10 - ' -0.2 

n 24+ 004 

„ „„+0.05 
U -' '-0.05 

11.9j*° 
73.3I3 ^ 



CMBall 

0226+ 0007 
u.uzzu_ 0006 

n in«+ - 007 
U.1U8_ 006 

1.042l° ;»I 
0.0921°:°" 
0.82l°it 

n 960 +00ie 

u.you_ 015 

o Q4+0.06 

0.75+ 003 



13.7Z 



-0.03 
7+0.1 
-0.2 

n 2^+ 03 

u - zd -0.03 

0.78+ 004 
11.6 
0.92 



0.04 
+2.3 
2.7 
+0.08 
-0.15 

72.9l£g 



CMBall+BIMA CMBall+BIMA+LSS 

0226+ 0006 
u.uzzo_ 0006 

n ins+ 007 
u.iu»_ 006 

1 n42+ 003 

l.U-*Z_ 003 

n nqi+o oi4 
u.uyi_ 014 

70+0 06 

u - ' 3 -0.09 

u.yoi_ 014 
o n4 +o.07 

0.76«» 3 

13.6±S;1 
0.24«» 3 
0.78^ 

11A+_H 
Q2+ 05 

u - 3z -0.06 

73Tlfj 



9 0226+ 0006 

u.uzzu_q 0006 

n 114+0.005 

i.042±S:ggi 

0.090lHl 2 2 

o.77t° ;2 

0.958lHl 5 5 

°- uo -0.05 

7 3+o.02 
13.7^ 



-0.03 
7+0.1 
-0.1 
n 97 +0.03 
u - z ' -0.02 
1+0.03 
-0.03 
+2.0 
2.4 
+0.05 



0.81 
11.5 
0.90 
70.7 



+2.1 




Fig. 1. — The de-correlated ACBAR band-powers for two alternate binnings. These two binnings are not independent, therefore only one 
set is shown with error bars, which correspond to 1-<t uncertainties calculated from the offset lognormal fits to the likelihood function. Both 
the overall features and the damping scale are in good agreement with predictions from a fiat, low baryon density ACDM Universe. The third 
acoustic peak (around I = 800) is clearly seen. Small scale primary CMB fluctuations are detected with high signal-to-noise ratio (> 4) up 
to I = 2000. The plotted model line is the best fit to the WMAP3 and ACBAR bandpowers. 



This factor R includes the full dependence of the calibration on RCW38's spectrum and the bandpasses of each exper- 
iment. The domi nant source of uncertainty is RCW38's spectrum; to be conservative, we double the estimates listed in 
iMasi et all ((2006). The model consists of two components: a power law term with a = 0.5 ± 0.2 and a dust term with 
Tdust = 22.4 ± 1.8 K. Only the relative amplitude of the two terms is important: A power i aw @3oaHz/Ad ustp eak = 867 ±400. 
We also include uncertainty in the laboratory measurement of each experiment's bandpass. The mean value and uncer- 
tainty in R is estimated using 100,000 realizations of the above parameters, and found to be 1.008 ± 0.021 (See Spectral 
Correction in Table lA"6l ) Given that our integration radius is larger then RCW38's size, the flux contribution of diffuse 
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Fig. 2. — Systematic tests performed on the ACBAR data. Top: Power spectrum produced from the difference of maps made with left- and 
right-going chopper sweeps. Middle: Power spectrum (diamonds) calculated from the difference of maps made from the first and second halves 
of the data for each field (CMB2, CMB5, CMB6, and CMB7), compared with Monte Carlo simulations (error bars). Bottom: Power spectrum 
calculated from the LMT sum maps (squares), compared with the joint power spectrum (pre-decorrelated, filled circles). The consistency 
between two sets of band-powers demonstrates that the residual chopper offsets are below noise. See text for more detail on this test. 



emission near RCW38 can be significant. The spectrum of this extended structure may be different from that of RCW38, 
in which case the calibration ratio would depend on the integration radius. We estimate this uncertainty from the observed 
variability of the calibration ratio with integration distance. 

The calibration value from the real map is normalized by the spectral correction for RCW38 and the signal-only transfer 
functions estimated for each experiment. The result of this analysis is that the temperature scal e for ACBAR's CMB 
fields in 2002 should be multiplied by 1.128 ±0.066 relative to the planet-based calibration given in lRunvan et a l. (2003). 
Table lA"6l tabulates the contributing factors and error budget. We now proceed to propagate this RCW38-based calibration 
to the CMB2 observation done in 2001. 



ACBAR 2001-2002 Cross Calibration We carry the 2002 RCW38 calibration into 2001 by comparing the 2001 
observations of the CMB2 field to the overlapping 2002 CMB4 field. The fields are reduced to the overlapping region 
and a power spectrum is calculated for each field. The bands are widened (A£ ~ 200) to avoid large noise correlations 
between the band-powers. Care was taken to insure the filtering of the two maps only occurs in the overlapped region. 
However, differences in scan patterns and array configurations between the two seasons cause differences in filtering. We 
assume that the band-powers from field a and field (3 are q/ (/ = a, (3). If the relative calibration factor between fields 
a and (3 is rj, we can find the value r]Q that maximizes the likelihood function: 



£(77) oc exp 



(Qa,i - VQ0,i) 2 



a 



D. 



(Al) 



for band i, which can be found by Monte Carlo analysis. A separate 
rjo is an unbiased estimator, and to calculate its uncertainty. 
Vo, 



The quantity <r 2 D . is the variance of (q a ,i — qp,i 
Monte Carlo simulation was carried out to confirm that 

We find the calibration factor to be CMB2/CMB4= 1.238 ± 0.067 {y/rjo, in units of temperature). Approximating the 
uncertainties as Gaussian, it implies the CMB2 temperature scale should be multiplied by 0.911 ± 0.074 relative to the 
scale used for the analysis of K04. 

The day-to-day relative calibration for the 2002 CMB fields is determined using the measured flux of RCW38. The 
procedures used are outlined in more detail K04. During some parts of 2002, RCW38 observations are only available 
for one of the two rows of 150 GHz channels. We derive the relative calibration between two rows of bolometers during 
these periods using the CMB power spectrum comparison method described earlier. We find the corrections to the 
BOOMERANG-based calibration factors are 1.031 ± 0.025, 0.935 ± 0.050, and 0.998 ± 0.042 for CMB 5, CMB6, CMB7, 
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Fig. 3.— The AC BAR band powers plotted with those from WMAP3 ijHinshaw et alJl2006ft and the 2003 flight of BOOMERANG 
ijJones et al.ll2006h . The three experiments show excellent agreement in the overlapped region. 

ACBAR and CBI excess 




Fig. 4. — ACBAR results on the high-£ anisotropics. Left: The ACBAR band-powers at I > 1000, plotted on a logarithmic scale with 
the latest CBI data taken at a frequency of 30 GHz. All the ACBAR bins at i > 2000 are lower than the CBI band-power measurement. 
Right: The likelihood distribution for the ratio of the "excess" power, observed by CBI at 30 GHz and ACBAR at 150 GHz. The excess 
for each experiment is defined as the difference of the measured band-powers and the model band-powers at t > 2000. The vertical dashed 
line represents the expected ratio (4.3) for the excess being due to the SZ effect. If the excess power seen in CBI is caused by non-standard 
primordial processes, the ratio will be unity (blackbody), indicated by the dotted line. We conclude that it is 4.5 times more likely that the 
excess seen by CBI and ACBAR is caused by the thermal SZ effect than a primordial source. In addition, because of the weak detection of 
excess power in ACBAR (1.2<r), it is about 3 times more likely that the excess is due to the SZ effect than radio source contamination of the 
lower frequency CBI data, assuming no contaminations from dusty proto-galaxies. 



respectively. We apply these corrections, and determine the overall calibration uncertainty to be 6% (in temperature 
units) based on the uncertainties associated with B03/ACBAR-2002 RCW38 cross calibration. 
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Fig. 5. — Basic parameter marginalized 1-dimensional likelihood distributions for the following data combinations; WMAP3 only (black, 
solid), ACBAR + WMAP3 (red, dashed), CMBall (green, long-dashed), and CMBall + LSS (blue, dash-dot). All runs include lensing. 
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Fig. 6. — Two-dimensional marginalized distribution for the erg and the spectral index n s . The contours are for the 68% and 95% confidence 
levels. The underlying (red) region corresponds to the WMAP3 only basic parameter set in both left and right panels. In the left panel the 
overplotted (yellow) region is for the combination of ACBAR + WMAP3. In the right panel we show the combination ACBAR + CMBall. 
The Harrison-Zel'dovich-Peebles, scale invariant solution with n s = 1 is only ruled out at the 2cr level. 




Fig. 7. — Two-dimensional marginalized distribution for the correlated pair Urun = dn s /dlnk and n s for the three data combinations. As 
in Fig. \E\ the contours are for the 68% and 95% confidence levels. The straight dashed lines show the scale invariant case. This illustrates 
that negative running is preferred at less than a 2-cr level. If a tensor component was added, the errors on running would become even larger. 
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Fig. 8. — Two-dimensional marginalized distribution for the conventional <rg parameter (left panels) and SZ template fit derived <r| z (right 
panel), plotted vs. the spectral index n s . The contours shown are for the 68% and 95% confidence levels. The data combinations used are 
shown in the insets. The horizontal lines show the la confidence region for the <rg parameter obtained from the 1-d marginalized posterior for 
the CMBall + BIMA + LSS combination with no SZ template fitting. The SZ derived values for the normalization of the matter fluctuations 
are higher than those found from the basic parameter run. 
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Error Budget for the RCW38 based ACBAR Calibration 
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Spectral Correction 
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Note. — The calibration of ACBAR through RCW38 has multiple factors and potential sources of error, tabulated here for 
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Fig. A9. — On the left-hand side is a map of RCW38 made by B03. The bandpasses for each experiment measured by Fourier transform 
spectroscopy are shown on the right. The range of possible spectra (ltr) for RCW38 is also shown, with each spectrum being normalized to 
0.50 at 150 GHz. (green and violet lines). 
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